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ABSTRACT 

We present high signal-to-noise measurements of the acoustic scale in the presence of nonlinear 
growth and redshift distortions using 320h~^ Gpc^ of cosmological PM simulations. Using simple fit- 
ting methods, we obtain robust measurements of the acoustic scale with scatter close to that predicted 
by the Fisher matrix. We detect and quantify the shift in the acoustic scale by analyzing the power 
spectrum: we detect at greater than 5a a decrease in the acoustic scale in the real-space matter power 
spectrum of 0.2% at z = 1.5, growing to 0.45% at 2 = 0.3. In redshift space, the shifts are about 25% 
larger: we detect a decrease of 0.25% of at z = 1.5 and 0.54% at z = 0.3. Despite the nonzero amounts, 
these shifts are highly predictable numerically, and hence removable within the standard ruler analysis 
of clustering data. Moreover, we show that a simple density-field reconstruction method substantially 
reduces the scatter and nonlinear shifts of the acoustic scale measurements: the shifts are reduced to 
less than 0.1% at z = 0.3 — 1.5, even in the presence of non-negligible shot noise. Finally, we show 
that the ratio of the cosmological distance to the sound horizon that would be inferred from these fits 
is robust to variations in the parameterization of the fitting method and reasonable differences in the 
template cosmology. 

Subject headings: distance scale — cosmological parameters — large-scale structure of universe — 
methods: N-body simulations — cosmology: theory 



1. INTRODUCTION 

Baryon acoustic oscillations (BAO) imprint a dis- 
tinctive feature in the large-scale structure of both 
the cosmic microwave background and the matter 
density field. The sound waves propagating through 
the primeval plasma of photons and baryons in 
the early Universe freeze out when the photons are 
freed from baryons at the epoch of recombination. 
This l eaves unique oscilla t ory features in the CMB 
(e.g.. iMiller et all 119991: Ide Bernardis et all l2000l: 



2003t iBlake fc Clazebrool [20031 iHu fc Haimal] [20031 : 
Seo fc Eisensteinll2003f ). Observing the BAO from large- 
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set by the distance that the sound waves have traveled 
before recombination, known as the sound horizon 
or acoustic scale. Since the sound horizon scale can 
be well determined by present and future CMB data, 
the BAO can serv e as an excellent sta n dard ruler 
(|Hu fc Whitd [1996[ : [Eisenstein et al.[ [1998[ : [Eisensteid 
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scale structure at different redshifts with the knowledge 
of their true physical scale allows us to geometrically 
measure the expansion history of the Universe and 
thereby identify the driving force behind t he observed 
accelerated expansion, i.e., d ark energy (e.g.. [Riess et al.l 
[19981 [Perlmutter et al.[[l999[ l. 

Using the BAO as a standard ruler for precision cos- 
mology requires that we understand all of the physical 
effects that could alter the acoustic scale during the evo- 
lution of the Universe. There are two important aspects 
of performance to consider for a standard ruler measure- 
ment with BAO: precision and accuracy. The precision of 
the standard ruler method crucially depends on the non- 
linear evolution of structure growth and observational 
effects, such as redshift distortions and galaxy bias. 
With time, these nonlinear effects increasingly degrade 
the contrast of the BAO in the matter power spectrum 
and correlation function, decreasing the signal-to-noise 
of the standard ruler method. Recently, significant 
progress has been made on modeling the nonlinear 
evolution of the amplitude of the BAO with analytic 
and/ or numerical meth o ds ([Meiksin. White, fc Peacockl 
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I2006t ISeo fc EisensteinI l2007f) . We can now make 
reasonable predictions for the degradation of the BAO 
due to nonhnear growth and redshift distortions. 

The accuracy of the standard ruler method depends 
on the calibration of the physical scale of the BAO, i.e., 
the sound horizon. Future BAO surveys aim to push 
the standard ruler method to near the cosmic variance 
limit. This will require a calibration of the sound horizon 
scale to about 0.1% accuracy. However, nonlinear evolu- 
tion not only degrades the precision of the standard ruler 
method by decreasing the contrast of the BAO, but also 
may alter [i.e., shift] the observed BAO scale in the large- 
scale structure at low redshift. Such a shift, relative to 
the linear acoustic scale derived from the CMB, will de- 
grade the accuracy of the ruler. Failure to appropriately 
account for such a shift will result in biased measure- 
ments of the cosmological distances and therefore bias in 
the inferred d ark energy parame t ers. 

Previously, iSeo fc EisensteinI (|2005[ ) had no detec- 
tion of a shift (or bias) in the acoustic scale at the 
~ 1% level in the nonlinear real-s pace pow er spec- 
trum from N-body data (also see iHuff et al. 2007 ; 
Mai 120071 lEisenstein et "all 120071 I Seo fc Eisenstcin. ,2007: 
Angulo et all l2008l : ISanchez et a l. 2008). Meanwhile, 
Smith et all (|2008f ) predicts more than a percent level 
shift in the acoustic scale, when defining the observed 
scale as the local maxima of the peak in the correlatio n 
function (also see lSmith et al.ll2007t [Cuzik et al.ll2007D . 
However, the effect of the broadband shape, as well as 
the effect of the nonlinear smoothing of the BAO, on the 
shift of the acoustic scale can be largely marginalized 
over with appropriate parameterizations, which is prob- 
ably the reason why suc h a l arge shift is not dete cted in 
ISeo fc EisensteinI ()2005f l and lAngulo et all (|2008D . That 
is, the choice of acoustic scale estimator matters; one 
cannot quantify a nonlinear shift without defining the 
estimator and one would surely prefer to use one that is 
close to optimal. The latter goal recommends template 
fitting to a significant portion of the power spectrum or 
correlation function, rather than peak finding methods. 

Recent analytic and numerical work by 
ICrocce fc Scoccimarrol (|2008D makes use of renor- 
malization perturbation theory and predicts a nonzero, 
irreducible shift at the few tenths of a percent level 
in real space due to nonlinear mode coupling. This is 
referred to as a "mode-coupling shif t"; a similar contex t 
is introduced as "physical shifts" in lSmith et al.l (|2008l ). 
The exact level of this mode-coupling shift would also 
depend on estimators: for example, the effects of the 
mode-coupling shifts that are different on different nodes 
may average out to some degree in the case of a global 
fit over several acoustic peaks. The good estimators 
will be ones that marginalize over most of the effects 
from the broadband shape and minimize contributions 
from the mode-coupling shifts. Any residual shifts can 
be modeled through numerical computation, such as a 
study of a large volume of A^-body data. 

In this paper, we show that most of the acoustic 
scale information can be extracted via a template fit- 
ting method. In order to test and confirm such a low 
level mode-coupling shift in redshift space as well as in 
real space, we require a very large volume of cosmolog- 
ical simulation. We use 32Qh~^ Gpc^ of simulations to 
test the precision and the sub-percent level accuracy of 



the standard ruler method with BAO. In this paper, we 
consider both nonlinear growth and redshift distortions. 
The effect of galaxy bias will be studied in our subsequent 
p aper (Siegel et al . , in p reparation) . 

'Eisenstc in et al.l ()2007l ) show that a portion of the non- 
linear degradation of the BAO can be undone with a sim- 
ple reconstruction sch eme based on the ZeP dovich ap- 
proximation (also see ISeo fc EisensteinI I2007D . We test 
how well we can undo the nonlinear shift as well as the 
erasure of the BAO, even given the presence of shot noise. 

In § [21 we describe our cosmological N-body simula- 
tions and methods of analysis to measure the acous- 
tic scales from the simulations. In § |3l we present the 
resulting shifts and errors on the measurements of the 
acoustic scale when accounting for the nonlinear growth. 
In § |4j we proceed to redshift space, and present the ef- 
fects of redshift distortions. In § jSj we show that not 
only the nonlinear erasure but also the nonlinear shift on 
the acoustic scale can largely be undone via reconstruc- 
tion. In § |6l we compare our results of signal-to-noise 
with the Fisher matrix calculations. In § [71 we test how 
our results depend on the small errors on the template 
cosmology used in our fitting. Finally, in § |8l we summa- 
rize the major results obtained in this paper and point 
towards future directions for probing dark energy using 
BAO as a standard ruler. 

2. SIMULATIONS AND THE METHODS OF 
ANALYSIS 

We have chosen to investigate a specific cosmology of 
the ACDM family: = 0.25, = 0.75, f^b/i^ = 

0.0224, h = 0.7, n = 0.97 and cts = 0.8. In order to model 
the non-linear evolution of structure, we make use of N- 
body simulations. Since our focus is on large scales and 
we need a lot of volume, we have elected to run 40 real- 
izations of a 2 /i~^Gpc box using a parallel particle- mesh 
(PM) code. Each simulation employs 1024'^ particles on 
a 2048^^ force mesh, for a mesh resolution of lft.~^Mpc 
and a particle mass of 5 x IQ^^^ h~^MQ. The acoustic 
scale is thus resolved by ~ 100 grid cells in each dimen- 
sion and the non-linear scale by ^ 10 cells. The initial 
condition s were generated using the Zcl'dovich approx- 
imation (jZel'dovichl Il970f ) starting from a regular grid 
at z = 50. The linear theory power spectrum for the 
initial conditions was computed by evolution of the cou- 
pled Einstein, flui d and Boltzmann equa tio ns using the 
code described in IWhite fc Scott I (|1996t ). ISeliak et all 
(120031) find that this code agrees well with CMBfast 
(jSeliak fc Zaldarriagalll996l ). We used constant steps in 
In a with steps of either 5% or 10%. The larger time 
steps resulted in a loss of power at small scales, but we 
saw no dependence on time step in the acoustic scale 
information. 

Power spectra were computed at z = 1.5, 1.0, 0.7, and 
0.3 by assigning the mass to a grid using the cloud- in-cell 
interpolation and using wavenumber bins of width Afc = 
0.0047 h Mpc~^. In order to reduce the data volume that 
we need to store and manipulate, we save only 1% of the 
processed data for the later use (§[5]). In this paper, we 
use this 1% fraction only for our studies of reconstruction 
and shot noise. 

We perform a analysis to fit the spherically aver- 
aged power spectrum Pobs(fc) to template power spectra 
Pmik/a). We parameterize the observed power spectra 
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as 



Pobs(fc) = B{k)P^{k/a) + A{k), 



(1) 



2005 



2006 



2008 



where a, Bik), and A{k) are fitting parameters. Here 
a is a scale dilation parameter and represents the ra- 
tio of the true (or linear) acoustic scale to the measured 
scale. For example, a > 1 means that the measured 
BAO being shifted toward larger k relative to the lin- 
ear power spectrum. In the case we fail to account for 
any nonlinear shift on the acoustic scale in the standard 
ruler method, a represents the ratio of the mis-measured 
distance to the true distance. The term B{k) allows a 
scale-dependent nonlinear growth, and A[k) represents 
an anomalous power, i.e., additive terms from the non- 
linear growth and shot noise. By including both B[k) 
and A{k) with a large number of parameters, we mini- 
mize the contribution to the standard ruler method from 
the broadband shape of the power spectrum. We try var- 
ious parameterizations for B{k) and A{k), as described 
in §|S]and §|4l 

We calculate the model or template power spectrum 
Pm by modifying the BAO portion of the linear power 
spectrum with a nonlinear parameter Em to account 
for the degradation of the B AO due to non linear 
effects and red s hift d ist ortions (lEisenstein et al.l 
Tegmark et al.l 120061: ICrocce fc Scoccimarrd 
Eiscnstein et aLl 120071 : ICrocce &: Scoccimarrd 
Matsubara 200l): 

Pm{k) = [Fii„(fc) - Fsmooth(fc)] exp [-eY?j2\ 

+Psmoothik), (2) 

where Pun is the linear power spectrum and Psm ooth is 
the nowiggle form from lEisenstein fc Hul ()1998D . Any 
details of physics overlooked in the template power spec- 
trum will be, at least partly, absorbed into B{k) and 
A{k). In fact, we will show that our results do not 
vary for a wide range of Em, as the fitting formula can 
alter the amplitude of the BAO by trading power be- 
tween B{k)Pm and A{k), as B{k) and A{k) have enough 
flexibility. Fine-tuning the choice of Em in Pm does 
not affect our results, and therefore we do not include 
Em as a fitting parameter. Wc use a fitting range of 
0.02/1 Mpc-i < fc < 0.35/1 Mpc-i. 

As the true covariance matrix of the A^body simula- 
tions is unknown a priori, we use the variation between 
the simulations to assess the true scatter in a. We fit 
each simulation assuming the covariance matrix for P(fc) 
is that of a Gaussian random field, i.e., assuming inde- 
pendent band powers with variances determined by the 
number of independent modes in each. While this is not 
an optimal weighting of the data for the determination 
of a, the effects from non-Gaussianity in the density field 
will still be reflected in the scatter between best-fit a's 
from different simulations. 

We use different groupings to assess the scatter be- 
tween our simulations. Our primary approach is to mea- 
sure the scatter in a from jackknife sampling of the 40 
simulations. However, we have also used bootstrap sam- 
pling for 40 sets of 1 simulation. The results from the two 
approaches are highly consistent; we will quote results for 
the jackknife case. Note that the scatter in a, i.e., ctq, 
that we quote represents the scatter associated with the 
mean value of a for a total volume of 320/i~^ Gpc'^, not 
the scatter between samples. 
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Fig. 1. — Real-space (top) and rcdshift-space (bottom) power 
spectrum P(k) from our A'^-ljody simulations for z=0.3, 0.7, 1.0, 
and 1.5. The power spectra are di vided by Psmooth i the nowiggle 
form from lEisenstein fc Hul II1998I ') , scaled with D'^ . The dashed 
line is the linear power spectrum. The gray lines indicate the large- 
scale amplitude expected from linear theory. 



3. A REAL-SPACE ANALYSIS OF THE BAO 

3.1. Non-linear Growth of Structure 

Nonlinear growth arises when the density perturba- 
tion on a given scale reaches an amplitude of order 
unity; the evolution of perturbations of different wave 
modes becomes increasingly coupled with one another, 
causi ng a substantial d e parture f r om linear evolution 
'e._g,^ lJuszkiewiczl 119811: Ivishniad 119831: iMakino et al l 



iggj iJain fc Bertschingd 119941: rBharadwail ll996allbl: 



Meiksin fc Whitelll999l : IScoccimarro et al.lll999f) . Such a 
nonli near process erases the BAO of the power spectrum 
(e.g. iMeiksin. White, fc Peacockll 19991: ISeo fc EisensteinI 
l2005f ). thereby degrading the signal of the standard ruler 
method. It also increases small-scale power above the 
linear growth rate, thereby increasing the noise in the 
clustering statistics. This nonlinear process first appears 
on small scales and proceeds to larger scales with time. 
From the perspective of configuration space, the non- 
linear growth damps the BAO peak at ~ 100/i~^ Mpc 
because the large-scale bulk flows cause the differential 
motions of the dens ity pairs initially sepa rated by the 
sound horizon scale (jEisenstein et al.ll2007l ). In this sec- 
tion, we measure the effect of nonlinear structure growth 
on the BAO using the real-space power spectra from our 
A-body simulations. 

3.2. Fitting Results in Real Space 

In the upper panel of Figure [1] we illustrate the 
spherically-averaged real-space power spectra from our 
A-body simulations for each of our four redshifts (0.3, 
0.7, 1.0, and 1-5), div ided by the nowiggle form from 
lEisenstein fc Hul ()1998l ). As expected, the effect of non- 
linearity increases the small-scale power and increasingly 
degrades the BAO at lower redshifts. 
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Wc measure the shifts and the scatters of the shifts 
of the BAO relative to hnear theory with a analysis. 
We model the real-space nonlinear power spectrum by 
parameterizing B{k) and A{k) in equation ([T]) in two 
independent ways. We use polynomial fitting, which we 
refer to as "Poly?" , with 

"Poly 7" : B{k) =bo + bik + 62/fc^ 

^(fc) = ao + aifc + a2fc^ + ... + ayfc'^, (3) 

as well as fitting to Pade approximants. which we refer 
to as "Pade" , where 



"Pade": B{k)^bo 



1 + cik + c^k'^ + cr,k^ 



1 + C2fc + c^k"^ 
A{k) = Oq + aik + a2k^. 



(4) 



Here , 6i , and Ci are free fitting parameters for all i . In 
addition, we consider a case in which A{k) in equation 
([3]) is truncated at third order in k, which we refer to as 
"PolyS" . In constructing P,„ . we use reasonable choices 
of Sm that are based on the Zel'dovich approximations 
(|Eisenstein et al.ll2007[ ): 7.6/1"^ Mpc at z = 0.3 and scal- 
ing with the linear growth factor at other redshifts. We 
vary Sm near our fiducial values so as to show that the 
analysis is not sensitive to the choice of Em. 

We perform a jackknife analysis on our 40 simulations 
of P{k), using the two different fitting forms above. We 
detect non-zero shifts on the mean acoustic scale [i.e., 
a — 1] of a few tenths of a percent at all redshifts, as 
shown in Table [1] These shifts represent a decrease in 
the fitted acoustic scale [i.e., the fitted BAO being at 
larger k], as opposed to the linear scale, of 0.20%, 0.26%, 
0.33%, and 0.45% for respective redshifts of 1.5, 1.0, 0.7, 
and 0.3. The shifts increase at lower redshifts, as ex- 
pected from the evolution of nonlinear structure. Our 
results are statistically significant at the 5f7, 6a, 7cr, and 
8ct levels for respective redshifts of 1.5, 1.0, 0.7, and 0.3, 
for both fitting forms. The reduced values in the ta- 
ble are obtained from the best fit to the averaged power 
spectrum over 40 simulations and show the goodness of 
the fit in general. The reduced x^i X^/DOF. is close to 
unity even though we have used a Gaussian approxima- 
tion for the errors on the nonlinear power spectrum. 

Our results are consistent over a wide range of Em 
(AEm = zt2h~^ Mpc) used in Pm, as shown in the lower 
half of Table [TJ the difference in a is well within ctq, and 
the difference in ctq is well within the expected Gaus- 
sian sample variance of x 40 « 11% for 40 simula- 
tions. For example, varying Em between 6.Qh^^ Mpc and 
9.0h~^ Mpc at z = 0.3 causes negligible changes in a and 
tTa (less than 10^'^ in a). Reasonable variations at other 
redshifts are similarly negligible. That is, due to the large 
number of parameters allowed in our parameterization, 
our fitting is able to tolerate a mild discrepancy between 
Em used in P,„ and the actual power spectrum, without 
a substantial degradation of the fit. 

While the "Pade" model gives a slightly better fit than 
the "Poly7" model, based on the smaller x^/DOF val- 
ues, we note that both models give the same acoustic 
peak scale and very similar errors. We tried varying 
the order of the fitting formulae and found that our 
results did not change. Of course, one cannot make 
B{k) and A{k) arbitrarily flexible or else these terms 
would be able to mimic the BAO. With a fitting range 




0.2 

k (h Mpc-') 

Fig. 2. — Real-space power spectra at z = 0.3 and 1.5 divided 
by the initial linear power spectrum: wc divide Pres = {Pni ~ 
A{k)) / B{k) by Pun. The oscillatory feature represents the amount 
of degradation of the BAO. Red error bars arc Gaussian errors 
centered at Pros/Piin- Black solid lines are for the template power 
spectra Pm{k / a) / Pii^{k) constructed by using — 7.%h~^ Mpc 
and 4.5h~^ Mpc at 2 = 0.3 and 1.5, respectively. Blue solid lines 
represent the nowiggle form, i.e., Psuioothi^ / Oi) / P\\n{k) ■ 



of 0.02/1 Mpc-1 < fc < 0.35/1 Mpc-\ the flexibility we 
allowed in both our fitting formula avoids this problem. 
As shown in Table [1] decreasing the polynomial orders 
of A{k) to the third order, i.e., "Poly3", makes a negli- 
gible change to the best-fit a and a a- However, using a 
model in which B{k) is constant produced notably worse 
fits and made the fit more sensitive to the choice of Em 
used in Pm ■ To be conservative, we recommend a second 
order Bik) and higher than fourth order in A{k) in real 
space. 

We find that the bootstrap analysis is highly consistent 
with the jackknife estimate of the scatter in a. Both 
agree to within 5% in ctq and 0.0001 in a. 

Figure [2] shows the real-space average power spectrum 
aX z ~ 1.5 and 0.3, and, as a comparison, the best-fit 
models to them. Here, we first subtract the additive 
part A{k) from the nonlinear power spectra and then di- 
vide the residual Pa\ — A{k) by B{k) and by the input 
linear power spectrum, so as to manifest the degrada- 
tion in the BAO feature. That is, the broad-band differ- 
ences, relative to the initial power spectrum, have been 
removed and only the difference in the BAO structure is 
illustrated: any oscillations in the figure represents the 
nonlinear erasure of the BAO. We note that the BAO 
feature in Figure [2] depends on the value of Em we as- 
sume for the template Pm- The results at z = 0.7 and 1 
are similar to the redshifts shown. Error bars are drawn 
assuming Gaussian errors, which will be an underestima- 
tion on small scales, where there should be a contribu- 
tion from the trispectrum. We also ignore the nonlinear 
covariance between the different data points in the fig- 
ure. The excellent agreement between the best-fit power 
spectrum and the A'^body data demonstrates that our 
template model has a suitable amount of flexibility to 
match the data well. 
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Fig. 3. — Left and Middle panels: the growth of nonlinear shift in q — 1 with redshift. The data points show the mean of a as a function 
of redshift; the errors show the scatter. The points are correlated. In the left panel, we overplot curves proportional to D^; in the middle 
panel, we overplot curves proportional to D^. Right panel: the a — 1 at z = 0.3 and z = 1.5 for 40 simulations. One can see that a given 
simulation produces a that are highly correlated between redshifts. 



To conclude, we find shifts of 0.2% to 0.45% in the 
acoustic scale due to nonlinear growth. The shifts in- 
crease at lower redshifts. Our recovery of the acoustic 
scale from the template fitting method is robust for dif- 
ferent fitting forms, a wide range of Sm, and the choice 
of jackknifc or bootstrap resampling. 

3.3. The Growth of a — 1 with Redshift 

We next investigate the redshift dependence of the shift 
a. The left and the middle panel of Figure [3] shows the 
growth of a — 1 with redshift with the error bars ctq, 
for 32Qh~^ Gpc'^. The dotted lines are curves of lin- 
ear growth with different proportionality constants in 
the left panel, and are curves scaling quadratically with 
the growth factor in the middle panel. From the figure, 
a — 1 is growing faster than the growth factor, D, but 
slower than . We find that a — 1 scales as I?™ where 
m = 1.66^q'22, when we perform a fitting to the mean 
a with a covariance matrix generated by the variations 
in the bootstrap samples. The right panel shows that 
values of a — 1 for each subsample are highly correlated 
between different redshifts, as expected. A model of a — 1 
scaling as D^, w hich one might expect from the per- 
turbation theory (jCrocce fc Scoccimarrol [20081 ) . is only 
1.3a away from the best fit; we cannot reject this model. 
Meanwhile, a — 1 scaling as D is disfavored at about 3a. 
However, a model with leading order D and higher-order 
terms fits well and cannot be excluded. 

We next compare our results with the analytic predic- 
tions. While we detect the sub -percent leve l nonh near 
mode-coupling shift predicted bv lSmith et al.l (I2008D and 
Crocce fc Scoccimarrol ()2008f ) (also see iNishimichi et al.l 
2007( ). the redshift-dependence of a — 1 we derive is 
weaker than the node-by-node sh ifts in the BAO feature 
in pow er spectrum predicted by ICrocce fc Scoccimarrol 
(l2008f). where they si mulate the fitting method used in 
iPercival et all (|2007bl ) ^ Note, however, that we are us- 
ing a global fit. where the overall shift will be a result of 
different nodes being weighed differently depending on 
redshift, resulting in an overall redshift dependence that 
is not the same as the redshift dependence of individual 
nodes. Indeed, when they consider a global fit, the red- 
sh ift dependence of the shift seems attenuated [Eq. (27) 
of ICrocce fc Scoccimarro (2008,)]. In addition, other de- 
tails of our fitting method are different from the cases 

^ Their a corresponds to the inverse of a defined in this paper. 
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Fig. 4. — Redshift-spacc power spectra at z = 0.3 and 1.5 divided 
by the initial linear power spectrum: we divide P^^s = (^nl ^ 
A{k))/[B{k)F{^g] by Pu^- The oscillatory feature represents the 
amount of degradation of the BAO. Red error bars are Gaussian 
errors centered at Prea/Pun- Black solid lines: Pm,{k / a) / Pii^{k) 
constructed by using Em = 9.0/i~^ Mpc and 6.0/i~^ Mpc at z = 
0.3 and 1.5, respectively. Blue solid lines represent the nowiggle 
form, i.e., Psmoothik / a) / Pun{k) ■ 

they considered. Taking these differences into account, 
we consider that these analytic works and our A'^body 
results are largely in agreement. 

4. A REDSHIFT-SPACE ANALYSIS OF THE BAO 

4.1. Redshift- space Distortions 

The observed power spectrum is subject to rcdshift- 
space distortions due to the peculiar velocity of the ob- 
served sources. The large-scale bulk motions [i.e., outflow 
from undcrdcnsc regions and infall toward overdcnse re- 
gions] along the line of sight increase the apparent density 
contrast in this direction and induce an angle-dependent 
amplifi cation of the real-space powe r spectrum on large 
scales (iKaise^ 119871: iHamiltonl fl998t IScoccimarroll2004l: 
iPapai fc Szapudil 120081 etc). The apparent displace- 
ments due to such redshift distortions degrade the BAO 
further in redshift space relative to real space (e.g. , 
iMeiksin. White, fc Peacockll999l : [Se"o fc Eisensteinll2005l : 
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lEisenstein et al.l l2007f ). Meanwhile, virial motions be- 
tween and within halos cause an apparent suppression 
of the small-scale power along the li ne of sight, which is 
called the finger-of-God effect (FoG) (|de Lapparent et al.l 
[l98l . 

While the redshift distortions arc an inevitable element 
of the study of cosmological distance error estimation, 
a direct analysis of the redshift-space power spec- 
trum has been challenging due to the lack of a reliable 
model for fitting redshift distortions and the anisotropic 
nature of the redshift-space power spectrum (cf . , 
Wagner et all l2007l: lA^eulo et al.l[200l lOkumura et all 



2008 HPadmanabhan fc Whitdl2008l) . Recent theoretical 



and n umerical studies (jEisenstein et al.ll2007l : iMatsubaral 
l2008f ). however, provide quantitative predictions for the 
BAO degradation due to nonlinearity and redshift distor- 
tions and therefore for the resulting distance errors (e.g., 
ISeo fc Eisensteinll2007[ ). Based on these studies, we can 
improve the description of the nonlinear effect on BAO 
in the fitting model and check the distance errors derived 
from the redshift-space power spectra. 

In this section, we extend the analysis to the 
spherically-averaged redshift-space power spectrum to 
find the effect of the redshift distortions on the acous- 
tic scale measurements. We do not include large-angle 
effects in the redshift distortions but merely use the flat- 
sky approximation. Due to a computer hardware prob- 
lem, some of the redshift-space data were lost, so that 
our redshift-space analysis uses a volume of 216/i~^ Gpc^ 
(i.e., 27 simulations). 

4.2. Fitting Formulae 

The lower panel of Figure [1] shows the spherically aver- 
aged redshift-space power spectra from our simulations. 
In order to account for the FoG effect in a analysis, 
we add a simple multiplicative function F^g to our pa- 
rameterization. We try two ways of incorporating Ff^g 
into our parameterization: 



^'In": P{k) = B{k)P,^{k/a)F,,g+A{k), 



(5) 



and 



"Out": P{k) = [B{k)Prn{k/a) + A{k)] x Ffog 

+ ei. (6) 

where A{k) — ao + aik + a2k^. Note that the highest 
order of polynomials for A{k) is only k^, less than that 
of "Poly 7" in real space. We lower the number of param- 
eters in A{k) in redshift space, as there are enough addi- 
tional parameters to account for the broad-band shape. 
The parameter ei is added to permit a shot noise term 
in the "Out" model . 

We again try both polynomial and Fade approximants 
for B{k): 



■Tolf: B(fc) = 6o + 6ifc + 62fc^ 



(7) 



and 



We": i?(fc)=6o ^+r'+r-'''tr'' . (8) 



1+C2fc + C4fc2 



We try two forms for Ftog: 



"Rat" : Ffog = 1/ [l + (fcdi)''^] (9) 



with priors of di > O.lh ^ Mpc and ^2 > 0, and 

"Exp": f^fog = exp[-(fcdi)'^^] (10) 

with priors of di > Oh~^ Mpc and d2 > 0. We refer 
to, for example, the combination of equations ([7]) and 
pO|) in the form of equation ([6]) as "Poly-Exp-Out" . In 
addition, we also try "Poly7" (eq. [3] ) to fit the redshift- 
space power spectrum. We will mainly quote the results 
from using "Poly-Exp-Out" . 

We increase the values of T,^ used in P,„ relative to 
the choices in real space, so as to account for the addi- 
tional degradation of the BAO due to redshift distortions. 
Note that we are approximating the anisotropic nonlin- 
ear effect on the BAO in redshift space with an isotropic 
parameter Sm. 

4.3. Fitting Results in Redshift Space 

Table[2]lists the errors and the mean of the best-fit a of 
jackknife samples in redshift space. We also present the 
values of CTq rescaled to a survey volume of 320h~^ Gpc^ 
from the actual results for 216h~'^ Gpc^ so as to provide 
an easy comparison to Table [1] 

We again detect a non-zero, sub-percent shift of the 
mean values of a in redshift space: 0.25% of shift of the 
acoustic scale at z = 1.5 by 5a, 0.33% at z = 1 by 5a, 
0.41% at z = 0.7 by 6ct, and 0.54% at z = 0.3 by 6a, 
where the signal-to-noise ratios are calculated based on 
the original cTq for 216h~^ Gpc'^. The shift on a has 
systematically increased relative to the real-space values 
at all redshifts. We find that a increases by ^ 25% at all 
redshifts; we do not detect a redshift-dependence in this 
increase. As expected, values of CTq also have increased 
relative to the real-space values at all redshifts: about 
25% at z = 0.3 and 0.7, 18% at z = 1, and 11% at 
z = 1.5. 

We find our results to be consistent for both a and CTq , 
regardless of the various parameterizations. As shown in 
Tabic [2] for z = 0.3, the differences in a are well within 
(Jq and the differences in CTq, are well within the sample 
variance, 1/^2 x 27 « 14%, in most of the cases. Inter- 
estingly, we find that even without correcting for i^tog, 
such as the real-space fitting formula "Poly 7" (eq [3]), 
yields results consistent with ones with Ffog correction, 
so long as enough free fitting parameters are allowed. 
Again, the results arc not very sensitive to the fiducial 
values of S,n used to calculate Pm- We therefore choose 
equation ^ combined with the polynomial B{k) and the 
exponential form of Ffog (i.e., eq [10]), "Poly-Exp-Out", 
as our convention to fit the redshift-space power spectra 
for the rest of paper. 

Figure 2] shows the redshift-space power spectrum at 
z = 1.5 and 0.3 and the best-fit models to them. Again, 
we first subtract the additive part A{k) from the nonlin- 
ear power spectra, but this time we divide the residual 
Pa\ — A{k) by B{k)F{og before dividing by the input lin- 
ear power spectrum, so as to manifest the degradation 
in the BAO feature. The figure illustrates the excellent 
agreement between the best-fit power spectrum and the 
A-body data. 

5. RECONSTRUCTING THE BAO 

In this section, we reconstruct the baryon acous- 
tic signature obscured by nonlinear growth and red- 
shift distortions using the simple scheme presented in 
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Fig. 5. — Rcal-spacc P-ni(k) of the 1%-samplc after reconstruction (red lines), divided by the nowiggle form from lEisenstein fc Hul II1998I ). 
in comparison to the nonUnear power spectra before reconstruction (blue lines). The gray lines are for the large-scale amplitude expected 
from linear theory. The dashed lines are for linear power spectrum. The effect of reconstruction is most apparent at z = 0.3 and is smallest 
at 2 = 1.5. 
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Fig. 6. — Real-space Pn\{k) of the 1%-sample after reconstruc- 
tion (red error bars): we divide Pres = (^^nl ~ ^i^)) / B(k) by 
Plin- The oscillatory feature represents the amount of degrada- 
tion of the BAO. Red error bars are Gaussian errors centered at 
Pros/Piin- Black solid lines are for the template power spectra 
Pm{k / ct) / Pi[ii{k) constructed by using hsted in Table [3] Blue 
solid lines: Psmooth(fc/")/-Plin(fc)- 

lEisenstein et al.l ()2007f ). We wish to investigate whether 
such reconstruction impro ves the signal-to- noise ratio 
on the measurement of a (ISeo fc EisensteinI l2007f ) and 
whether it reduces the shift of the acoustic scale. 

For the analysis with reconstruction, we use the phase- 
space data of a randomly selected 1% of all the particles. 
We will refer to this sample as "1%-sample", in compar- 
ison to the previous "100%" sample. For the 1% sam- 
ple, power spectra are computed for wavenumber bins of 
width Ak = 0.001ft. Mpc~^ and we use a fitting range of 
0.02/1 Mpc"^ < k < 0.5/iMpc"^ We utilize a volume 
of 320h~^ Gpc^ both in real space and in redshift space. 
We note that by using the 1% sample, we perform the 
reconstruction scheme in the presence of non-negligible 



shot noise. ISeo fc EisensteinI ()2007f ) quantify the ratio 
of shot noise to sample variance with nPo.2, where n is 
the number density and Po.2 is the real-space power at 
k = 0.2h Mpc~^. For the analysis here, nPo.2 runs from 
0.7 at z = 1.5 to 2.4 at z = 0.3. 

We apply the linear perturbation theory continuity 
equation to the nonlinear density fields of the 1%-sample 
and predict the linear theory motions. We first take the 
nonlinear density fields at z=0.3, 0.7, 1.0, and 1.5 and 
Fourier-transform them. We apply a 10h~^ Mpc Gaus- 
sian filter to S{k) in order to smooth gravity on small 
scales. We then compute the linear theory motions q us- 
ing V ■ 9 = Finally we displace the real particles and 
smoothly distributed reference particles by the resulting 
q, and find the final dens ity field from the differ ence of 
the two density fields. As lEisenstein et all ()2007[ ) noted, 
this differencing method moves the measured densities 
back to their initial locations: linear structure growth is 
mostly preserved and only the nonlinear growth effect is 
reduced as a result. 

5.1. Reconstruction in Real Space 

Figure [5] shows the reconstructed power spectra of the 
1% sample, divided by Psmooth, at various redshifts. The 
benefits of reconstruction are most obvious at z = 0.3 
due to the low shot noise and the strong nonlinearity to 
correct for at this redshift. We fit the power spectra of 
these reconstructed density fields using the same models 
as before (eq. [3] and [4]). We find consistent results for 
both models in general. Because the upper limit of the 
fitting range was chosen to be fc = 0.5h Mpc""'^ in the 1% 
case, we find that the "Poly 7" model with more fiexible 
A{k) fits better to the reconstructed real-space power. 
We decrease for the template P,„ approximately by 
half relative to Sm values used before reconstruction. 

Table[3]shows the resulting a—1 and CTq, for the "Poly 7" 
model. Because the reconstruction analysis uses the 1% 
sample, the results in Table [3] are not immediately com- 
parable to those in Table [l] which used the 100% sample 
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and hence had negUgiblc shot noise. Wc include results 
for the 1% sample without reconstruction in Table [3] to 
provide the proper comparison. 

We find that applying reconstruction decreases Ua by 
almost a factor of 2 at z = 0.3. The improvement is 
incrementally less at higher redshifts. This is to be ex- 
pected for several reasons. First, the density field is more 
linear at high redshift. As most of the acoustic scale in- 
formation \s at k < 0.2ft, Mpc~^, the non-linear degrada- 
tion of a a is smaller at high redshift. There is therefore 
less for the reconstruction algorithm to fix. Second, the 
intrinsic clustering is weaker at high redshift such that 
the shot noise is more important. This contaminates 
the measurement of the density field and hence the esti- 
mation of the reconstruction displacements, causing the 
reconstruction to be less accurate. Third, the increased 
shot noise means that the high harmonics are not being 
well measured, regardless of whether or not they are cor- 
rectly reconstructed. This means that the results with 
and without reconstruction will not differ much. For ex- 
ample, the comparison of Tables [T] and [3] (summarized in 
the "N-body" column of Table [5]) shows that the intro- 
duction of shot noise degrades by over a factor of 2 at 
z = 1.5, even without reconstruction. Finally, the force 
resolution of our PM simulations {Ih^^ Mpc mesh) may 
not be adequate to recover the BAO with full fidelity 
when our reconstructed Sm is expected to be as small as 
2 — 2>h~^ Mpc, as it is at z > 1. If our small-scale struc- 
ture is more fluffy, this will appear as an extra source of 

In addition to the improvement in the scatter of a. Ta- 
ble [3] also shows that reconstruction decreases the shifts 
in a itself. The shift in a decreases to 0.07% at z = 0.3 
and to below 0.03% at higher redshifts. Note that while 
reconstruction docs not improve ctq at z — 1.5, it does 
substantially improve the shift in a. This likely indicates 
that the reconstruction method is correctly compensat- 
ing for the large-scale bulk flows but is leaving a small- 
scale jitter that washes out the high harmonics of the 
BAO, as discussed above. Figure [5] shows the best-fit 
models to the averaged power spectra of the 1%-sample 
after reconstruction. Increasing the Gaussian filter size 
to \Ah~^ Mpc makes no difference: less than 2% in CTq, 
and 10^*^ in a. 

5.2. Reconstruction in Redshift Space 

Next we apply the reconstruction techniques to the 
redshift-space density field. The method is unchanged, 
save that we use the observed nonlinear redshift-space 
density field and apply the resulting displacement field 
to the redshift-space positions (ISisenstein et al. 2007). 
In this paper, we do not atte mpt to correct for the FoG 
effect of virialized halos (c.f.. lEisenstein et al.ll2007l ). 

Figure [7] shows the redshift-space power spectra after 
reconstruction, divided by the nowiggle form; the differ- 
ence before and after reconstruction is significant even in 
redshift space, not only at z = 0.3 but also at z = 0.7 
and 1.0. Table m quantifies the values of a and its scat- 
ter after applying the reconstruction method in redshift 
space. We find that the scatter is substantially re- 
duced to a level comparable to that found with the re- 
constructed real-space case. That is, the reconstruction 
method seems more effective in redshift space. The BAO 
is more smeared in redshift space, providing more oppor- 



tunity for reconstruction to improve CTq. In addition, the 
higher amplitude of the power spectra in redshift space 
reduces the effects of shot noise and may permit better 
reconstruction. 

As in real space, the shifts in the acoustic scale, a, 
are markedly improved, decreasing to below 0.1% at all 
redshifts. Figure [5] shows the resulting best-fit to the 
averaged power spectra, compared to the data. Again, 
we find a good agreement between the two. 

6. COMPARING THE N-BODY RESULTS TO THE 
FISHER MATRIX PREDICTION 

We next compare ctq from the A'-body data with the 
analytic distance error estimates from the Fisher matrix 
formalism. We calculate the Fisher matrix estimates 
by us ing a fitting formula based on ISeo fc Eisensteir] 
(|2007| ). This formula uses a parameter Eni to quantify 
the amount of non- linear degradation of the BAO. We 
estimate a reasonable value of Sni at each redshift by 
using the analyt i c Zel' dovich approximations (eq. [9] of 
lEisenstein et all (|2007l) ). 

Table [5] compares ctq from the N-body data with the 
Fisher matrix estimates. We start with the Fisher ma- 
trix errors corresponding to the 100% samples, for which 
shot noise is negligible; we use nPo.2 — 100 for all red- 
shifts. We do not include the nonlin ear growth effect 
into rt Pn.2, following the convention in lSeo fc EisensteinI 
()2007f ). While the Fisher matrix errors are in reasonable 
agreement with the N-body data, the Fisher matrix cal- 
culation systematically overestimates ctq by ~ 25%, 17%, 
15%, and 5% at z = 0.3, 0.7, 1.0, and 1.5, respectively, 
as shown in Table [5] Such an overestimation is con- 
trary to our intuition that the A-body data would reflect 
extra nonlinear effects, such as an effect on variance of 
the power spectrum, that the Fisher matrix does not in- 
corporate. As we allow many free parameters in our fit- 
ting forms, the contribution to the standard ruler method 
from the broad-band shape of the power spectrum is neg- 
ligible. However, because CTq is estimated from only 40 
simulations, the error on the scatter is expected to be 
1/V2 X 40 ~ 11%, assuming Gaussianity. While the dis- 
crepancy becomes larger than the expected sample vari- 
ance at z < 1, it is possible that the N-body results are 
lower than the Fisher matrix predictions simply because 
of sample variance. Because a values in a simulation are 
largely correlated between redshifts, it would not be sur- 
prising that the residuals in (Tq could all have the same 
sign. 

For the 1%-samples (Table [S]), shot noise is not negli- 
gible. We derive nPo,2 from the nominal shot noise, i.e., 
the inverse of the number density of particles. We find 
that the Fisher matrix errors are bigger than those from 
the simulations by 18% and 23% at z = 0.3 and 0.7, 
respectively. At z = 1.0 and 1.5, the agreement is bet- 
ter. The degradation in CTq from the simulation due to 
the addition of shot noise appears well predicted by the 
Fisher matrix formula. We note that the shot noise and 
its effect on a is not correlated between redshifts, leading 
to a little more scatter around the redshift trends. 

We next compare ctq, from the A'^body data with 
the Fisher matrix errors in redshift space. We 
use the 2-dimen s ional [2-D] approximation based on 
ISeo fc Eisensteir] (|2007[ ) and contract the 2-D matrix to 
derive the 1-D distance error, as the transverse and the 
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Fig. 7. — Redshift-spacc P{k) of the 1%-samplc after reconstruction (red lines), divided by the nowiggle form, in comparison to the 
redshift-space nonlinear power spectra before reconstruction (blue lines). The gray lines are for the large-scale amplitude expected from 
linear theory. The dashed lines are for linear power spectrum. The effect of reconstruction is most apparent at z = 0.3 and is smallest at 
z = 1.5 

A''-body data in real and redshift space are somewhat 
smaller than the Fisher matrix errors with Sni based 
on the Zel'dovich approximations. The measured er- 
rors from the iV-body data correspond to reductions of 
Sni by as much as lh~^ Mpc relative to the Zel'dovich 
prediction. It is possible that this systematic discrep- 
ancy is merely due to sample variance, which is ex- 
pected to be correlated between redshifts. However, 
there could be other reasons for the difference. It is pos- 
sible that the Zel'dovich approxi mation slightly overesti- 
mate s the non-linear degradation (jCrocce fc Scoccimarrol 
l2006l Padmanabhan & White, in preparation) or that 
the approximations behind the Fisher matrix calculation 
are slightly pessimistic. Previo us comparisons to high- 
resolution N-body simulations ()Seo fc EisensteinI l2007f) 
found agreement to better than 13%, comparable to the 
precision tested here. In principle, we could try to use the 
amplitude of the wiggles in P{k) to estimate Sni, but our 
template fitting parameter Em is not a precise proxy for 
Sni because the marginalization parameters allow good 
fits for a range of Sm- 

In conclusion, we are encouraged by the good agree- 
ment between the N-body and Fisher matrix results. The 
fact that the CTq, from the N-body fits are not much worse 
than those from the Fisher matrix indicates that our es- 
timator is extracting most of the acoustic information. 

7. EFFECTS OF VARIATIONS IN THE TEMPLATE 
COSMOLOGY 

When the true acoustic scale is known, the observa- 
tions of the BAO allow us to measure the cosmological 
distance scale, in particular the angular diameter dis- 
tance Da{z) and the Hubble parameter H{z). The true 
acoustic scale is determined from measurements of the 
matter and baryon density, e.g., from CMB anisotropy 
data. In this paper, we do not marginalize over uncer- 
tainties in the cosmological densities but instead use a 
fixed fiducial template. Using a fiducial template power 
spectrum assumes that a true acoustic scale as well as the 
shape of the linear power spectrum is known a priori. 
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Fig. 8. — Redshift-spacc power spectra of the 1%-sample after 
reconstruction (red error bars), in comparison to the best- fit models 
(black lines) and nowiggle form (blue lines): we divide Pres = 

(Pnl - A(fc))/B(fc)/Ffog by Plin. 



line-of-sight distance scales are the same in this simu- 
lated Universe. We use an ellipsoidal smoothing kernel 
with Exy (= E„i) and = [1 -I- rim(-z)°'^]Sxy and in- 
clude the angle dependence in nPa,2- The FoG effect 
is not included in nPo 2, as justified in lSeo fc EisensteinI 

(gool). 

From Table O the CTq, values from the 100%-samples 
deviate from the Fisher matrix errors approximately by 
1.5 times the sample variance. Again, they are systemat- 
ically lower than the Fisher matrix errors. Note that we 
expect the sample variance of x 27 w 14% rather 

than 1/V2 X 40 in this case. For the 1%-samples, the dis- 
crepancy reaches up to 27% at z = 0.3, while the discrep- 
ancy is below the sample variance of 1/^/2 x 40 w 11% 
at high redshift. 

Overall we find that the distance errors from the 
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Here we investigate how the results depend on small 
deviations in the template power spectrum. We consider 
variations in f^m^^, ^th'^, and the spectral tilt n of 2%, 
chosen to be somewhat larger than the expected uncer- 
tainties from Planck CMB anisotropy data. When one 
uses a different cosmological template compared to the 
true simulation input, one of course measures a shift in 
a. However, this shift in a is very close to that predicted 
by the ratio of the sound horizons in the two cosmologies, 
where the sound horizo n is given by equation (6) in 
lEisenstein fc Hul (|1998[ ). We find that a/rs is recovered 
to no worse than ±0.02% within the above cosmological 
variations, as shown in Table [6l 

This implies that the measured acoustic scale is consis- 
tently determined, independent of reasonable variations 
in the template cosmology. Applying this standard ruler 
to cosmological distance inferences, one will consistently 
measure DA{z)/rs and H{z)rs despite variations in the 
template used. Of course, when one tests a given cosmol- 
ogy against a suite of cosmological measurements, one 
must account for the change in in using the BAO con- 
straints, just as one must self-consistently alter any other 
predictions for cosmological observables. 

8. CONCLUSIONS 

Baryon acoustic oscillations in large-scale structure can 
serve as an excellent standard ruler to probe the accel- 
eration history of the Universe. Such a standard ruler 
method is based on the premise that the characteristic 
scale of the feature, the sound horizon scale, is well mea- 
sured from CMB and remains fixed throughout the evo- 
lution of the Universe. 

In this paper, we have measured the shift in the 
acoustic scale with high signal-to-noise ratio using 
320ft-~^ Gpc^ volume of the PM simulations. We allow a 
large number of parameters in our fitting formula so that 
the standard ruler method is dominated by the BAO. 

We find nonlinear growth decreases the measured 
acoustic scale by less than 1%: we find shifts of 0.20% 
at z = 1.5, 0.26% at z = 1, 0.33% at z = 0.7, and 0.45% 
at z = 0.3, all determined at 5 to Scr. In detail, we 

find the redshift dependence of a — 1 scales as D ■°°-o.22^ 
where D is the linear growth function. The perturba- 
tion theory prediction of a scaling is consistent at 1.3 
a. We consider different fitting formula, polynomial or- 
ders, and resampling methods and find highly consistent 
results between various choices. 

Our detection of a sub-percent-level shift in the 
acoustic scale agrees with the analytic nodc-by- 
node predictions of nonlinea r mode-coupling shift by 
ICrocce fc Scoccimarrd ()2008| ) within the same order of 
magnitude. The agreement is better when we compare 
our results with the shifts they predict assuming a global 
fit. 

Redshift distortions increase the shift by ^ 25% rela- 
tive to real space. Using 21Qh^^ Gpc'^ of simulations, we 
detect a shift of 0.25% at z 1.5, 0.33% at z = 1, 0.41% 
at z = 0.7, 0.54% at z = 0.3 at the level of 5 to 6tT. 

Although our results confirm the previous analytic pre- 
diction that nonlinear growth and redshift distortions in- 
duce a small mode-coupling shift on the acoustic scale, 
we find that the shift can be reduced to less than 0.1% by 
undoing the nonlinear growth using a simple reconstruc- 



tion scheme (lEisenstein et al.ll2007l ). even in the presence 
of non-negligible shot noise. 

We compare our results from TV-body simulations with 
the analytic Fis h er m atrix error forecasts based on 
ISeo fc EisensteinI ()2007l ). We find the standard ruler 
method from the simulations and Fisher formalism fore- 
casts arc in agreement within a modest level of sample 
variance. In detail, we find that CTq. from the simulations 
is systematically better, by no more than 27%, than the 
Fisher matrix calculations. While this trend could have 
an underlying physical cause, it is consistent with sample 
variance, correlated between redshifts, from our limited 
number of simulations. 

We show that small variations in the template power 
spectrum around the true fiducial cosmology do not af- 
fect our results oi a/rg- Therefore, the cosmological dis- 
tance ratios DA{z)/rs and H[z)rs inferred from the mea- 
surements of the acoustic scale are robust to reasonable 
errors in the template cosmology. 

We argue that it is correct and appropriate to use the 
linear theory power spectra as a template for extract- 
ing the acoustic scale. The acoustic oscillations are not 
a perfect harmonic sequence nor perfectly exponent ially 
damp ed. For example, because of Silk damping (jSilkl 
|1968( ). small-scale perturbations in the baryons decou- 
ple from the CMB earlier than large-scale perturbations, 
implying that the sound horizon is not a single quan- 
tity. However, the correct waveform is computed in great 
detail in the Boltzmann codes. Using an approximate 
waveform will simply risk a loss of precision and accu- 
racy. The exact relation of the acoustic signature and 
acoustic scale in the CMB to that in the low-redshift 
large-scale structure depends upon the details of recom- 
bination, but this is exquisitely constrained by CMB ob- 
servations. Once one has accepted the assumptions in 
interpreting the CMB anisotropics to measure the densi- 
ties needed to set the sound horizon, it is straight-forward 
to use the predicted waveform for the low-redshift power 
spectrum. 

Finally, we emphasize that the nonlinear shifts on the 
acoustic scale are predictable numerically, as shown here 
for each estimator. With such a prediction, the shifts 
can be properly modeled and therefore can be removed 
in the standard ruler analysis of actual clustering data, 
even without reconstruction; the results can be cross- 
checked for consistency with the results after reconstruc- 
tion, when available. 

We have successfully computed the effects of nonlin- 
ear growth and redshift distortions on the acoustic scale 
using a large volume of simulated Universe. The galaxy 
power spectrum is expected to exhibit somewhat larger 
shifts than the matter power spectrum, as galaxies would 
sample more nonlinear regions. We will investigate this 
in a separate paper (Siegel et al., in preparation). En- 
couraged by the successful reduction of shift by a recon- 
struction technique applied to the mass distribution with 
a nonzero shot noise, we will also investigate the use of 
reconstruction of BAO in the galaxy power spectrum in 
future work. 
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TABLE 1 

The mean and the error of q in real space 



z 


Model 




a — i(7oj 




X /IJOl 


0.3 


Poly7 


?.6 


0.45 


0.055 


1.09 


U.o 


r^aae 


1 .0 


U.4o 


U.UO I 


1 nn 
i.UU 


0.3 


PolyS 


7.6 


0.43 


0.056 


1.24 


0.7 


PolyT 


6.3 


0.33 


0.046 


1.10 


n 1 




D.O 


U.oo 


U.U4D 


U.y i 


1.0 


Poly? 


5.5 


0.2? 


0.041 


1.13 


1.0 


Fade 


5.5 


0.26 


0.042 


0.99 


1.5 


Poly? 


4.5 


0.20 


0.03? 


1.1? 


1.5 


Pade 


4.5 


0.20 


0.03? 


1.04 


0.3 


Pade 


6.0 


0.45 


0.05? 


0.98 


0.3 


Pade 


?.o 


0.45 


0.05? 


0.99 


0.3 


Pade 


10.0 


0.45 


0.05? 


1.03 



Note. — Fitting range: 0.02/i Mpc"^ < fc < 0.35/i Mpc"^ (a 
total of ?0 data points). Reduced x^, xV^OF, where "DOF" 
is the number of degrees of freedom, is calculated from the best 
fit to the averaged power spectra over 320/i~'^ Gpc'^. 



TABLE 2 

The mean and the error of a in redshift space. 



z 


Model 


Em OL 


- 1(%) 




xVdof 


0.3 


Poly-Exp-Out 


9.0 


0.54 


0.070 (0.085) 


1.26 




Poly-Exp-In 


9.0 


0.55 


0.070 (0.086) 


1.25 




Poly? 


9.0 


0.54 


0.068 (0.083) 


1.24 




Poly-Rat-In 


9.0 


0.50 


0.084 (0.103) 


1.46 




Poly-Rat-Out 


9.0 


0.53 


0.069 (0.084) 


1.23 




Pade- Rat-In 


9.0 


0.55 


0.077 (0.093) 


1.28 


0.7 


Poly-Exp-Out 


8.0 


0.41 


0.05? (0.069) 


1.12 


1.0 


Poly-Exp-Out 


7.0 


0.33 


0.049 (0.060) 


1.10 


1.5 


Poly-Exp-Out 


6.0 


0.25 


0.041 (0.050) 


1.11 


0.3 


Poly-Exp-Out 


8.0 


0.54 


0.071 (0.086) 


1.268 


0.3 


Poly-Exp-Out 


9.0 


0.54 


0.070 (0.085) 


1.256 


0.3 


Poly-Exp-Out 


10.0 


0.54 


0.070 (0.085) 


1.259 


0.3 


Poly-Exp-Out 


11.0 


0.5? 


0.070 (0.085) 


1.413 


Note. 


— Fitting range 


: 0.02/1 Mpc" 


-1 < fc 


< 0.35/1 Mpc"\ 


Values of 



Ga from the /V-body (the fifth column) are the errors for 40 simulations 
(320/i~^ Gpc^) rescaled from the errors for 27 simulations (216/i~^ Gpc^). 
The actual errors for 27 simulations are listed within the parentheses. Re- 
duced is calculated from the best fit to the averaged power spectra over 
216/i"^ Gpc^ 
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TABLE 3 

The effect of density-field reconstruction in real space. 



z Model Before Reconstruction After Reconstruction 







Em 


a - 1(%) 




xVdof 




Q - 1(%) 




xVdof 


0.3 


Poly7 


7.6 


0.40 


0.071 


1.09 


3.0 


0.07 


0.044 


1.15 


0.7 


Poly7 


6.3 


0.27 


0.062 


1.07 


2.0 


0.03 


0.049 


1.19 


1.0 


Poly7 


5.5 


0.22 


0.075 


1.00 


2.0 


0.00 


0.052 


1.13 


1.5 


Poly7 


4.5 


0.19 


0.079 


0.99 


2.0 


0.03 


0.080 


1.04 



Note. — Fitting range: 0.02/i Mpc"^ <k< 0.5/i Mpc"^ (a total of 481 data points). We use the 1% 
sample. 



TABLE 4 

The effect of density-field reconstruction in redshift space 



z Model Before Reconstruction After Reconstruction 







Em a 


- 1(%) 




xVdof 


Em Ot 


- 1(%) 




xVdof 


0.3 


Poly-Exp-Out 


9.0 


0.48 


0.079 


1.31 


4.0 


0.07 


0.046 


1.20 


0.7 


Poly-Exp-Out 


8.0 


0.30 


0.081 


1.14 


4.0 


0.03 


0.047 


1.21 


1.0 


Poly-Exp-Out 


7.0 


0.31 


0.093 


1.19 


3.0 


0.01 


0.056 


1.28 


1.5 


Poly-Exp-Out 


6.0 


0.27 


0.079 


1.08 


3.0 


0.04 


0.070 


1.19 


Note. 


— Fitting range: 


0.02/1 Mpc" 


^ < fc < 0.5/1 Mpc" 


We use 


the 1% sampl 


e. 







TABLE 5 

Fisher Matrix Estimates in Comparison to Simulation Results 



TV-body data Fisher matrix 











rr (%) 


E 1 

^nl 


71 Pn n 




Real space 


0.3 


100% 


7.6 


0.055 


7.6 


100 


0.068 




0.7 


100% 


6.3 


0.046 


6.3 


100 


0.054 




1.0 


100% 


5.5 


0.041 


5.5 


100 


0.047 




1.5 


100% 


4.5 


0.037 


4.5 


100 


0.039 




0.3 


1% 


7.6 


0.071 


7.6 


2.36 


0.084 




0.7 


1% 


6.3 


0.062 


6.3 


1.47 


0.076 




1.0 


1% 


5.5 


0.075 


5.5 


1.10 


0.075 




1.5 


1% 


4.5 


0.079 


4.5 


0.72 


0.078 




z 


Sample 




a„(%) 




nPo.2(M = 0) 




Redshift space 


0.3 


100% 


9.0 


0.070 (0.085) 


7.6/12.1 


100 


0.085 




0.7 


100% 


9.0 


0.058 (0.071) 


6.3/11.0 


100 


0.069 




1.0 


100% 


7.0 


0.049 (0.060) 


5.5/10.1 


100 


0.060 




1.5 


100% 


6.0 


0.041 (0.050) 


4.5/8.6 


100 


0.049 




0.3 


1% 


9.0 


0.079 


7.6/12.1 


2.36 


0.10 




0.7 


1% 


8.0 


0.081 


6.3/11.0 


1.47 


0.089 




1.0 


1% 


7.0 


0.093 


5.5/10.1 


1.10 


0.084 




1.5 


1% 


6.0 


0.079 


4.5/8.6 


0.72 


0.081 



Note. — Values of Em in the fourth column represents the nonlinear smoothing used for the template 
P(k) in the analysis of the iV-body data, and Eni and Exy/Ez in the sixth column are derived from 
the Zel'dovich approximations and represents the nonlinear degradation of the BAO assumed in the 
Fisher matrix calculations. For the 100% samples in redshift space, we list both the actual CTq from the 
27 simulations (inside the parentheses) and the rescaled for 40 simulations. 



TABLE 6 

The effect of a template P,n with an incorrect cosmology 





z 


Cosmology Alteration 


Ratio of a/rs 


Real space 


0.3 


1.02 VL-n^h^ 


1.0002 




0.3 


0.98 Q.^h^ 


0.9998 




0.3 


1.02 fift/i^ 


0.9999 




0.3 


Us = 0.95 


1.0000 




1.5 


1.02 ^mh^ 


1.0002 


Redshift space 


0.3 


1.02 VL-n^h^ 


1.0001 




1.5 


1.02 ^mh^ 


1.0002 



Note. — The ratio of a/rs, where is the sound horizon, in the 
fourth column compares the values of ajrs from an incorrect cosmol- 
ogy template with the one from the fiducial cosmology. Our fiducial 
cosmology is ^^h^ = 0.1225, Q.i>h^ = 0.0224, and n = 0.97. Fitting 
range: 0.02/i Mpc"^ <k< 0.35/1 Mpc~\ We use the 100% sample. 



